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ABSTRACT 

One of the most interesting results from the OGLE-III study of the LMC Cepheids is the large 
number of objects that pulsate simultaneously in the first and second overtone (denoted 10/20). 
Double-mode Cepheids yield important constraint on stellar evolution models. We show that great 
majority of the LMC 10/20 Cepheids have masses M = 3.0±0.5Me . According to current stellar 
evolution calculations, these masses are lower than needed for the blue loop in the helium burning 
phase to reach the instability strip. On the other hand, we found most of these stars significantly 
overluminous if they are crossing the instability before helium ignition. A possible solution of this 
discrepancy is to allow for a large overshooting from the convective core in the main sequence phase. 
We also discuss origin of double-mode pulsation. At the short period range we find two types of 
resonances that are conducive to this form of pulsation. However, at longer periods, it has a different 
(non-resonant) origin. 
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1. Introduction 

Double-mode (called also Beat) Cepheids are important objects for stellar physics. 
The diagnostic tool they provide are Petersen diagrams, where shorter to longer pe- 
riod ratio is plotted against log of the longer period. For nearly twenty years after 
they were introduced by Petersen (1973), the discrepancy between observational 
and calculated period ratio constituted a great challenge to stellar evolution the- 
ory. The major revision of stellar opacity calculation in the late 1980ties has been 
inspired by Simon's (1982) suggestion based mainly on this discrepancy. 

Only a small fraction of Cepheids develop double-mode pulsation. According 
to Soszynski et al. (2008b, hereafter S2008), it occurs only in less than 8 percent 
of the LMC Cepheids. It is difficult to develop this form of pulsation not only 
in nature but also in nonlinear modeling. Kollath et al. (1998) (see also Szabo, 
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KoUath and Buchler 2004) managed to obtain models exhibiting sustained double- 
mode pulsation but, as Smolec and Moskalik (2008b) argued, their treatment of 
turbulence was incorrect. 

So far attention has been mostly focused on Cepheids with excited fundamen- 
tal (F) and first overtone (10) modes. Readers are referred, in particular, to re- 
cent papers by Buchler and Szabo (2007) and Buchler (2008) on application of the 
F/10 Cepheid data as probes of stellar and stellar system metallicity. Surprisingly, 
data on the more frequent 10/20 double-mode pulsators did not attracted much 
attention. In the LMC there are more than three times as many 10/20 than F/10 
Cepheids (S2008). The present paper is devoted almost solely to the former type. 

2. OGLE-III data for the double- and triple-overtone Cepheids in the LMC 

The Petersen diagram for the 203 LMC 10/20 Cepheids is shown in the upper 
panel of Fig. 1. This is an expanded section of Fig. 2 in S2008, where also the 
P\o/Pf — log-Pp relation was shown. The patterns of these two relations differ. 
The latter shows a simple monotonic decrease. In our Fig. 1 we see a barely marked 
increasing tendency up to logPio ~ — 0.1 , followed by a sharp decline. The lower 
panel shows the period-luminosity relation (PL) for the same objects. The relation, 
in which the Wesenheit index is used as a reddening-free measure of luminosity, is 
part of the general PL relation for the first overtone Cepheids in the LMC depicted 
in the bottom panel of Fig. 6 in S2008. 

In each panel, we put also five points corresponding to triple-overtone objects 
(Soszyriski et al. 2008a). The three at logPio ~ -0.25 are 10/20/30 pulsators. 
The triple mode nature of two of them was first established with OGLE-II data 
by Moskalik et al. (2004). The three periods were then used by Moskalik and 
Dziembowski (2005) to construct seismic models of the object. The two F/10/20 
objects are widely separated in period. 

Soszynski et al. (2008a) found also two Cepheids with excited first and third 
overtone, showing no traces of the second overtone. These 10/30 pulsators have 
logPio = —0.277 and -0.244. The Wi indices ai^e within the band of the PL relation 
for the 10/20 Cepheids. 

Nearly half (101) of the 10/20 Cepheids occur in the \ogP\o range [-0.25,- 
0.05], that is, the 0.56-0.89 d period range, where it is the most common form of 
Cepheid pulsation in the LMC. Comparing histograms in Fig. 2, we see that there 
is much more of 10 Cepheids but not in this period range. There is no such a 
preferred range for the F/10 double mode Cepheids. In OGLE-III data for the 
LMC the number ratio of 10/20 to F/10 Cepheids is 3.2. The corresponding ratio 
from the EROS-2 data is only 1.8 (Marquette et al. 2009). This large difference is 
likely due to lower amplitudes of the 10/20 Cepheids and the lower peak detection 
threshold in OGLE-III dat£0. 

'We thank Igor Soszynski for this explanation. 
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Figure 1: The P20/P10 period ratio (top) and Wesenheit index, 
Wi = I — 1.55{V — I), (bottom) as function of the first overtone period (days) for 
203 10/20 (solid triangles), two F/10/20 (open squares), and three 10/20/30 
(open circles) Cepheids in LMC. 
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Figure 2: Period distribution of the 10/20 and 10 Cepheids in the LMC according 
to S2008 data. Total numbers of the 10/20 and 10 objects are 203 and 1228, 
respectively 

3. Masses, metallicity an evolutionary status of the 10/20 Cepheids in LMC 

Distribution of Cepheid periods is governed by various factors, such as recent 
history of star formation, time spend in the instability strip (IS), Xis , and mode 
selection. All these factors have been discussed by Alcock et al. (1999) in the con- 
text of the MACHO data on the LMC Cepheids. The value of Xis strongly depends 
on star mass and evolutionary status. For stars with intermediate masses, the first 
crossing of the IS takes place in the post-main sequence phase before helium ig- 
nites and it is by two orders of magnitude faster than subsequent crossings, which 
occur in the core helium burning phase (e.g. Alibert et al. 1999). Naturally, most 
of Cepheids are expected to be in the latter phase and have masses close to the 
minimum value, Mmis > for the blue loop of the evolutionary track in the HR dia- 
gram to reach up to the IS. In the OGLE-III sample containing over three thousands 
objects, we expect to find a significant number of stars during fast crossing of the 
IS. Furthermore, there are problems with explaining short period Cepheids as core 
helium burning objects (AUbert et al. 1999; Alcock et al. 1999). 

Estimate of Mmis is a subtle matter. The value depends on description of 
macroscopic element mixing, which is still uncertain. To illustrate the situation, we 
plot in the upper panel of Fig. 3 the evolutionary tracks downloaded from BASTI 
library (Pietriferini et al. 2004) for metalUcity parameter Z = 0.008, which is re- 
garded representative for young stars in the LMC. The tracks were calculated in 
two versions, adopting ttov = and 0.2 in the expression for the overshooting dis- 
tance X, = //pttov, where Hp is the pressure distance scale at the core boundary. 
The IS boundaries shown in Fig. 3 refer to simultaneous instability of the first and 
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Figure 3: The top panel shows evolutionary tracks from Pietriferini et al. (2004, 
BASTI files) calculated neglecting (solid lines) and including (dashed lines) over- 
shooting from convective core. The range of simultaneous instability of the first 
and second overtones (red dotted lines) is from our calculations. In the bottom 
panel the corresponding (only at M = 3 and 4Mq ) P-L relations are compared with 
the OGLE-III data for the 10/20 Cepheids. The two segments corresponding to 
the second and third crossing of the IS at M = 4Mq nearly coincide in this panel. 
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second overtones and they are shifted to higher temperatures by about 0.015 in 
log Teff relative to those for the F mode. The lower panel compares the PL relations 
corresponding to two versions of the tracks at stellar masses 3 and 4Mq with the 
OGLE-111 data for the 10/20 Cepheids. All pulsation properties were calculated 
by us (see below) for deep envelope models at the surface parameters taken from 
tracks and adopting the same input physics. Here and in the rest of this paper, the 
Wesenheit index for the models is calculated adopting 18.5 for distance modulus 
to LMC, which is close to the mean value from recent determinations by various 
methods (Shaefer 2008). The uncertainty of 0.05 of this value translates in less than 
0. IMq in inferred masses of the objects and has no bearing on our conclusions. The 
bolometric corrections and color indices are taken from Kurucz (2004). 

The Mmis = 4.5M0 value derived from the models calculated with ttov = 0.2 
is similar to that determined by Alibert et al. (1999) and Girardi et al. (2000). 
However, with ttov = the extent of the loop is not monotonic function of mass. 
Instead, it reaches maximum at M 4Mq and starts increasing again at M > 5Mq . 
Such a non-monotonic behavior is seen also at Z = 0.01 but not at Z = 0.004, and 
we do not know what is its origin. What matters for us here is that at Z correspond- 
ing to young LMC stars, with BASTI evolutionary tracks, we get helium burning 
Cepheids of shortest periods. Still, as we may see in the bottom panel of Fig. 3, 
by far most of the 10/20 Cepheids in LMC have shorter periods. Are they then 
objects which cross the instability strip for the first time? There are also problems 
with this hypothesis. At the short periods, the first crossing models with masses 
M < 3Mq and moderate overshooting seem consistent with data. However, this 
is not true for majority of the objects. Note that models of the same mass lay on 
nearly the same line in the logP — Wj plane, regardless of crossing number and the 
overshooting distance. We may see that most of the stars have masses closer to 3 
than 4Mq and they are overluminous. One may consider higher ttov value but the 
question is why such objects are so numerous though lis would be over 5 times less 
than for the short period objects and over 30 times less than for the 4Mq objects in 
helium burning phase. 

It would be easy to identify the IS crossing if we could measure rate of pe- 
riod change arising from evolutionary changes. Furthermore, such measurements 
would provide valuable constraints on stellar models. Unfortunately, the fast pe- 
riod changes, which are most often found in overtone Cepheids, have clearly non- 
evolutionary origin (Poleski 2008). 

The use of the second overtone periods allowed us to constrain parameters of 
the 10/20 Cepheids. In our study, we partially rehed on the old Warsaw codes 
calculating stellar models and their linear nonadiabatic oscillation properties (see 
e.g. Pamyatnykh 1999). However, these codes are not suitable for determination 
of the IS boundaries because convection is treated in a crude way. Therefore for 
this aim, we used a new code developed by Smolec and MoskaUk (2008a), which 
calculates radial mode frequencies and growth rates for unfitted envelope models 
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with a more sophisticated treatment of convection, wiiich allows to determine the 
range of the instability in a more or less reliable way. If the envelope is deep 
enough, the frequencies are found very close to those calculated using the old code 
for complete models with the same surface parameters. Our evolutionary models 
in the first crossing phase of the IS were used only as the reference. We varied 
luminosity relative to the evolutionary models, calculated with no overshooting and 
looked for the increment AL needed to explain periods of the 10/20 Cepheids. As 
the standard in our calculation we adopted the Asplund et al. (2002) heavy element 
mixture (hereafter AGS) and the OP opacity data (Seaton 2005). Convection is 
treated in the way described recently by Smolec & Moskalik (2008a) and as the 
standard, we adopt ttMcr = 1 -5 . We will briefly discuss consequences of alternative 
choices. There are more parameters describing time-dependent convection. Here, 
we use their values adopted by Baranowski et al. (2009, set B in Table 3) and they 
will not be discussed in the present paper. 

The main results of our survey are shown in Figs. 4 and 5, where calculated 
relations are superimposed on the data points. Let us first focus on the PL relations 
shown in the bottom panels. The effect of metallicity decrease is similar to that 
of luminosity increase, that is, the segments are moved upward along the constant 
mass fine. With only few exceptions, the points are contained between the 2.5 and 
3.5Mq Unes. We may see also that beyond the short period range, a significant 
excess of luminosity over the model values is required to match the data. The 
period ratio is indeed a sensitive probe of the excess. 

Most of the data points in the Petersen diagram fall into the range hmited by the 
segments corresponding to M = 2.5 and 3Mq , and luminosity excess AlogL = 0.4 
at Z = 0.006. Somewhat higher values of Z and/or M are needed to explain data 
at longest periods. Note in Fig. 5 that at Z = 0.008 the segment corresponding to 
M = 3.5Mq and AlogL = 0.4 runs through the mid of the data band for the longest 
periods. In Fig. 4 (Z = 0.006) such segment is not shown because it is outside the 
data range. In general, the models with Z = 0.006 are in better agreement with data 
but require larger A log L . Considerable spread of points at the shortest periods may 
perhaps reflect spread in the metallicity. 

The Pio/Pp period ratios are quite sensitive to Z. Using these data, Buchler 
(2008) derived ranges for Cepheid metallicities in various galaxies. For stars in the 
LMC he found the Z-ranges (0.0024, 0.0093) and (0.0062, 0.0124) with the AGS 
and GN93 (Grevesse and Noels 1993) mixtures, respectively. His calculations were 
done with the OPAL opacity data while in our survey we used the AGS mixture and 
OP data. Thus, it was important to check how use of different mixtures and opac- 
ities affects our conclusions. We found that the Pio/Pio dependence on P\o is 
much less sensitive to the choice of mixture than to opacity data. At low frequen- 
cies, where the effects are the largest the upward shifts of the P20/P10 ratio caused 
by use of the GN93 mixture is below 0.0005, while the downward shift caused 
by use of the OPAL data is about 0.002, which is similar to that caused by using 
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Figure 4: Calculated Petersen diagram (top) and PL relations (bottom) for selected 
models calculated with the metallicity parameter Z = 0.006 compared with the 
OGLE-III data for the 10/20 Cepheids. The segments correspond to the simultane- 
ous instability of the 10 an 20 modes for models of specified mass and luminosity. 
The red, blue, and magenta lines refer to 2.5, 3.0, and 3.5 Mq models, respec- 
tively. The solid lines correspond to evolutionary models calculated with ttov = . 
Dotted and dashed lines correspond to models with artificial luminosity increments 
AlogL = 0.2 and 0.4. Data are shown with solid triangles. 




Figure 5: The same as Fig. 4 but for models calculated with the metallicity param- 
eter Z= 0.008. 
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Z = 0.008 instead of 0.006. Thus, we do not see contradiction between Buchler's 
and ours assessment of Cepheid metallicity in the LMC. 

The choice of ttMLX has only a minor effect on P20/P10 ratio at specified P\o . 
The use of ttMLT = 1-8 instead of 1.5 causes the 0.0005 upward shift and a similar 
size downward shift is obtained with Kmlt = 10. Our conclusion regarding lumi- 
nosity excess holds for different choices of the opacity data, heavy element and the 
mixing length. However, because of the significant difference of the period ratios 
calculated with the two opacity data, we do not put much weight to our assessment 
of the Z values. 

Post-main sequence models with moderate overshooting and masses between 
2.5 and 3.5Mq explain data for the objects with logPio ^ —0.2. To this category 
belong the seismic models of the two 10/20/30 Cepheids constructed by Moska- 
lik and Dziembowski (2005). Certainly also the third such pulsator, which have 
has similar characteristics (see Fig. 1) may be interpreted as a standard post-main 
sequence object. The same appUes to the two F/10/20 pulsators. Note in that fig- 
ure that the star at logPio ~ 0.025 has the period ratio significantly higher than 
the 10/20 stars. However, to explain the observed characteristics of the great ma- 
jority of the 10/20 Cepheids in the LMC, we require a luminosity excess larger 
than that caused by overshooting within normally considered range, if indeed these 
objects are crossing the IS for the first time. From our code at M = 3Mq , we 
find AlogL = 0.55aov. The corresponding increment is by some 50% higher at 
M = 3.5Mq . From plots in Fig. 3, we may find that BASTI tracks yield somewhat 
higher increment resulting from overshooting and also that the difference in logL 
between the first and second crossing is about 0.25 at M = 4Mq . 

The large number of objects in the M = 2.5 and 3M0 mass range with luminos- 
ity excess AlogL = 0.4 is difficult to explain. If they are crossing the IS before he- 
Uum ignition then the excessive overshooting may mimic rotation induced mixing. 
Then, the problem to explain is why they are more numerous than the short period 
objects though, with larger luminosity, they spend less time in the IS. Alternative 
possibility is that the long-period 10/20 Cepheids represent a special population 
of heUum burning objects with relatively high fractional mass of the hydrogen-free 
core. 

4. Resonances and origin of double-mode pulsation 

The origin of double-mode pulsation is still lacking satisfactory explanation 
(Smolec and Moskalik, 2008b). Two ways in which the pulsation instabiUty may 
lead to such a terminal state were first discussed by Dziembowski and Kovacs 
(1984). One is a non-resonant way acting through the feed-back effect of pulsa- 
tion on mean radial structure. The alternative way involves resonant excitation of 
stable modes. A large number of papers has been devoted to the 2:1 resonance 
between the F and 20 modes. The coupUng to stable mode prevents the F mode to 
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Figure 6: The resonances during first crossing of tlie instability strip for stars with 
masses 2.5, 3.0, and 3.5Mq (red, blue, and magenta lines, respectively). Solid 
and dashed lines refer to models calculated with ttov = and 0.2, respectively. 
Positions corresponding to the V50 = 2vio resonance are shown with the asterisk 
and those to the 2v2o = Vio +V30 resonance with the pentagons. 
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saturate driving and allows excitation of the 10, which is linearly unstable. The ef- 
fect was seen in the numerical models but only with unreaUstic parameters (Kovacs 
andBuchler 1988, Smolec 2008). The three-mode resonance 2vio =Vf-|-V2o was 
shown by Smolec and Moskalik (2007) to cause double-mode pulsation in models 
of P Cephei stars. Unfortunately, we do not know any star of this type with only 
two radial modes excited. 

Both types of resonances, but involving higher order-overtones, occur in certain 
range of our models. In Fig. 6, we marked positions of the resonances within (and 
somewhat outside) of the IS in our models. Here, only complete evolutionary mod- 
els calculated with ttov = and 0.2 were used. We may see that at logPio ^ —0-3 
the equality V50 = 2vio occurs within the IS for our models. In fact a signature of 
resonance centered at logPio ~ —0.5 in the Fourier coefficient for the 10 Cepheids 
has been reported in S2008. Inspecting plots in Fig. 6, we may see that this implies 
masses around 2.5Mq for models calculated with ttov = 0.2 and somewhat higher 
if ttov is less. 

The resonant coupling to the strongly damped the 50 mode may cause a signif- 
icant amplitude reduction of the 10 amplitude allowing excitation of the 20, which 
is normally disfavored. This effect may explain origin of the double-overtone pul- 
sation in the short-period objects. At longer periods, the resonance occurs outside 
of the IS. We do not have complete models for most of the long -period objects so 
unfitted envelope models are unreliable for calculation of high overtone frequencies 
and one could speculate that the resonance might return to the IS in the range of the 
decreasing P20/P10 ■ However, we think that this is unUkely because we found that 
for the 4M0 BASTI models, which reproduce periods at the bottom of this branch, 
we are far from resonance. 

The 2v2o = Vio + V30 resonance occurs within the IS only in a narrow interme- 
diate period range. Thus, also this resonance caimot explain origin of double-mode 
pulsation in most of the 10/20 Cepheids. Such form of pulsation must predom- 
inantly arise from non-resonant saturation of the driving effect. Only nonlinear 
modeling may tell us how does it happen. The histogram in Fig. 2 yields an im- 
portant hint for such a modehng. Selecting a model with the first overtone period 
near 0.7 d, we are very likely to see the instability developing to double-overtone 
pulsation. 

The explanation in terms a non-resonant saturation seems difficult for the two 
10/30 Cepheids (Soszyhski et al. 2008). In these cases, the 2:1 resonance between 
the second and the seventh overtone may play a role. The resonant coupUng to this 
strongly damped mode may totally prevent excitation of the second overtone. 

5. Conclusion and discussion 

We confronted the OGLE-HI data on 203 LMC Cepheids, in which first and 
second overtones are excited (the 10/20 Cepheids), with model calculations. The 
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data used in our paper include periods and the Weseinheit index, Wi , which com- 
bined with the LMC distance is used as the measure of the absolute magnitude. 

Unlike the F/IO, the 10/20 Cepheids cannot be regarded an exceptional form 
of stellar pulsation. Not only the latter form occurs three times more often among 
the LMC objects but, more importantly, in a certain period range around 0.7d, it 
is the most common form. The data on the 10/20 Cepheids present interesting 
challenges to stellar evolution and pulsation theory. 

The period-luminosity (W? — logPio) relation constrains masses of the objects. 
Nearly all of them fall in the 2.5 to S.SMq mass range. The Petersen diagram 
{Pio/Pio — logPio) yields constrains on metallicity and luminosity. We focused 
on the former, adopting for the metallicity parameters Z = 0.006 and 0.008. The 
inferred luminosities for the objects with the first overtone period shorter than about 
0.6 d are well explained with post-main sequence stellar models crossing the insta- 
bility strip for the first time, calculated assuming no or moderate overshooting. 

However, for objects with longer periods, which constitute majority of the sam- 
ple, we find significant luminosity excess. If these are objects crossing the insta- 
bility strip before heUum ignition then the excess could be due to an overlarge 
overshooting, perhaps connected with fast rotation. The difficulty of this explana- 
tion is the short crossing time of the instability strip while the number of the objects 
showing the luminosity excess is relatively high. If these are helium burning objects 
then the difficulty is the low inferred mass. Standard evolutionary track for stars 
with M < 3.5 and acceptable metallicity do not enter the instability strip in this 
evolutionary phase. We must postulate a non-standard evolution leading to objects 
with relatively massive heUum cores. This could be a large mass loss in the red 
giant phase or again a large overshooting in the main sequence phase. The former 
option would be plausible if the stars were in binary systems for which we do not 
have yet any evidence. The latter option remains to be verified by means of new 
stellar evolutionary calculations. The result is not clear. As a warning, we have to 
recall that in certain mass range, BASTI models calculated without overshooting 
entered the instability strip in helium burning phase while those with overshooting 
did not (see Fig. 3). 

The challenge to stellar pulsation theory is to identify conditions leading to 
excitation of two rather than one pulsation mode and to explain why the 10/20 
pulsators are more frequent than F/10 and also 20. We found two types of reso- 
nances, which may promote double-overtone pulsation but only in the case of the 
short-period objects. For the majority of stars, such pulsation may result only from 
non-resonant saturation. In a certain period range this must be the preferred form of 
terminal development of the instability. The question why it is so may be answered 
only with the help of nonhnear modehng of stellar pulsation. 
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